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ABSTRACT 

Following this astro-ph communicat ion, we release a set of 10 000 WMAP ILC simulations, produced 
as described by lEriksen et alJ (J2004) . We strongly encourage that an analysis of these simulations 
accompanies any scientific analysis of the observed ILC map, in order to assess the stability of the 
results. For examples of such analyses, we refer the reader to our original paper. In particular, Figure 
lc) and d) show the bias and uncertainty induced by the ILC method, and Figure 7 shows a sobering 
examples of a \ow-t analysis; the reconstructed quadrupole-octopole alignment is plotted against the 
true input alignment, and the correlation is underwhelming. It can hardly be emphasized too strongly 
- the heavily processed full-sky WMAP maps are not reliable for quantitative analysis, neither at 
large nor small scales, and they should only, at best, be used as supportive evidence. We hope that 
these simulations will reduce some of the confusion regarding these maps that currently exists in the 
community. 

Subject headings: cosmic microwave background — cosmology: observations — methods: numerical 
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The main data product of the first-year WMAP data 
release llBennett et al]l2003a|) was a set of ten sky maps 
covering frequencies between 23 and 94 GHz. These sky 
maps formed the basis of both cosmological and astro- 
physical analyses, and unprecedented constraints on both 
the CMB power spectrum and cosmological parameters 
were obtained. 

For public convenience, a set of template-corrected 
sky maps were also provided by the WMAP team on 
the LAMBDA 1 web-site. For most cosmological appli- 
cations, these constitute the data set of choice (cou- 
pled with a proper sky cut), since they have both well- 
understood noise properties and fairly low levels of resid- 
ual foreground contamination. 

The WMAP team als o released a full-sky Int ernal Lin- 
ear Combination map l)Bennett et al.ll2003b|) intended 
for visualization purposes only. Its complicated noise 
and foreground properties (and mixing of the two) make 
it most unreliable for quantitative analysis. However, 
since the Galactic plane is not readily visible, it is useful 
for public display. The WMAP team explicitly warns 
against using it for CMB analysis. 

For some applications, full-sky maps are highly desir- 
able due to algorithmic limitations. In particular, meth- 
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ods that operate in spherical harmonics space are greatly 
simplified if no sky cut is imposed in the analysis. For 
this reason, the ILC map, which by eye looks almost free 
of foreground residuals, has been extensively used for sci- 
entific purposes - despite the fact that there are strong 
(and difficult to quantify) residual foregrounds present in 
the map. 

This provided the motivation fo r our first analysis o f 
the ILC map, which is detailed bv lEriksen et alJ {2004). 
In this study we found that the ILC map is indeed highly 
contaminated by residual foregrounds, and in particular, 
that the low-£ components, which have received the most 
attention so far, are highly unstable under the ILC clean- 
ing operation. A sobering example is the analysis of re- 
constructed versus input quadrupole-octopole alignment, 
which is summarized in Figure 7 of the above paper: 
Clearly, the preferred directions of the \ow-£ components 
are more or less randomized by the residual foregrounds 
present in the ILC map. 

During that analysis, we also found that the weights 
used to form the official WMAP ILC sky map were sub- 
optimal: Lower variances could be found by weighing the 
frequency bands slightly differently. This was due to poor 
convergence in the non-linear search algorithm used by 
the WMAP team, and was thus purely an algorithmic 
effect. Rather than relying on non-linear searches, we 
therefore solved the problem using Lagrange multipliers. 
However, since neither our map nor the one produced by 
the WMAP team is intended for quantitative analysis, 
the difference did not matter. But the two maps should 
not be considered independent, and the fact that they 
agree to some extent should not be taken as evidence of 
low foreground contamination. 

Given all the attention these maps have received, we 
believe it is useful to make a set of ILC simulations 
publicly available. These simulations are produced by 
adding foreground templates and channel-specific noise 
to random CMB realizations, and running these "semi- 
realistic" simulations through the ILC algorithm. By 
comparing the input CMB realization to the recon- 
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structed realization with some statistic of choice, one 
gains some intuition of the errors introduced by the ILC 
method. However, we point out that the errors are likely 
to be underestimated, because of the constant spectral 
indices implied by the foreground templates, and also 
because the template coefficients are estimated from the 
high-latitude sky only, and this leads to underestimation 
of the low-latitude foreground levels. These simulations 
should therefore be interpreted as setting a lower limit 
on the residuals present in the observed ILC map. 

It is also worth noticing that in order to estimate the 
errors in the ILC map, it is not sufficient to simply add 
scaled foreground templates to CMB realizations. For 
the WMAP ILC map, individual weights are computed 
for 12 disjoint regions, and the interplay between noise 
and foregrounds within each of these regions makes it 
very hard to predict the full-sky results. End-to-end sim- 
ulations such as the ones provided here is should be stud- 
ied to understand the effects of the ILC method on the 
statistic of interest. 

A set of 10 000 ILC simulations may be freely down- 
loaded here: 



http : / / www .astro . uio . no / ~hke / 

/cmbdata/L!LC_sim_input.tar.gz 



http : / / www. astro . uio . no / ~hke / 

/cmbdata/L!LC_rec_input.tar.gz 

The first file contains the input CMB realizations, and 
the second file contains the reconstructed ILC maps. 



For storage reasons, these two files contain the spher- 
ical harmonics coefficients ai m (in standard HEALPix 
format) up to £ = 50 for each map, which hopefully 
should be sufficient for more low-£ analyses. The 10 001. 
realization in the reconstruction set contains the actual 
observed LILC map. Other formats or ^-ranges may be 
obtained by sending an email to HKE. 

The spherical harmonics coefficients have been multi- 
plied with the response function of a Gaussian beam of 
1° FWHM and a pixel window function of corresponding 
to a HEALPix resolution parameter N s a c = 512. Thus, 
if a sky map is desired, as opposed to spherical harmonics 
coefficients, the data as provided here must be convolved 
with a beam of at least 8° FWHM to have a properly 
bandwidth limited function at ^ max = 50. If the func- 
tions are projected on grids with lower resolution than 
Aside = 512, one also must convolve with the appropriate 
pixel window. 

We also point out that if one wants to study the low- 
£ components of the WMAP data, the rec onstruction 
method described by iBielewicz et alJ l)2004fl applied to 
the template-corrected sky maps is significantly more 
accurate than the ILC method. The uncertainties in- 
duced by the cut sky are considerably smaller than the 
foreground-induced uncertainties in the ILC maps (see, 
e.g., Figure 1 of Bielewicz et al. 2005). 

We hope that this set of simulations will prove useful 
for future analyses, and that they may help reducing the 
confusion caused by these "foreground-cleaned" full-sky 
maps. 
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